Algorithms for uniform particle initialization in domains with complex boundaries

Accurate mesh-free simulation of fluid flows involving complex boundaries requires that the be captured accurately in terms particles. In context incompressible/weakly-compressible flow, SPH method is more accurate when particle distribution uniform. Hence, for time flow presence boundaries, one must have both an boundary discretization as well a uniform particles to initialize simulation. This process obtaining initial called "particle packing". this paper, various packing algorithms present literature are implemented and compared. An improved SPH-based algorithm proposed which produces distributions solid domains two three dimensions. Some challenging geometries constructed demonstrate accuracy new algorithm. The implementation open source manuscript fully reproducible.